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Abstract 

The formfactors of i?i — > vr and Di ^ tt transitions, where Bi{Di) is the 1^ 
P-wave Q^u^^q meson state, are calculated in the framework of the light cone QCD. 
Furthermore these formfactors are compared with the pole dominance model predic- 
tion using the values of the strong coupling constants QbiB'tv and QdiD^tt-, and a good 
agreement between these two different descriptions is observed. 
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1 Introduction 



Understanding the formfactors of the hadronic currents is one of the main problems in 
particle physics. As is well known, the knowledge of the formfactors allows one the de- 
termination of the quark mixing parameters, in a unique way [0]. Moreover, it gives us a 
direct information about the dynamics of the hadronic processes at large distance, where 
perturbation theory of QCD is invalid. At present there exists no theoretical method for 
calculating the formfactors, starting from a fundamental QCD Lagrangian. Therefore, in 
estimating the formfactors, usually phenomenological or semiphenomenological methods 
are used. 

Along these lines, there exists two well known methods for calculating the hadronic form 
factors. First of these is the pole model description based on the vector dominance idea, 
that suggests a momentum dependence dominated by the nearest pole, and the other one 
is the QCD sum rules method. The pole model approach does not have its roots on the 
basic principles of the theory and is purely phenomenological. Generally, it is assumed that 
the vector dominance approximation is valid at zero recoil, that is at —>■ m^, where is 
the square of the momentum transfer, that depends on the formfactor, and m is the mass 
of the nearest vector meson. However, there are no reliable arguments in favour of the pole 
model, for it to be also valid at small p^, that are interesting for practical applications. The 
effective heavy quark theory predicts somewhat larger region of validity, characterized by 
(m^ —p^)/mQ ~ (1 GeV), where mg is the heavy quark mass (for more detail see |0, H and 
the references cited therein). It is shown in 0, ^] that the pole behaviour of the formfactors 
is consistent with the dependence at — > predicted by the sum rule. This result has 
been confirmed independently by the calculations carried within the framework of the light 
cone sum rules ^ j^. 

In contrary to the pole model, the QCD sum rules method is based on the first principles 
and on the fundamental QCD Lagrangian. In this work we employ an alternative version 
of the QCD sum rules, namely light cone sum rules method for the calculation of the 
formfactors. We note that using this approach, the formfactors of the 7rA°7* transition 
P], the pion formfactors at the intermediate momentum transfer 0, semileptonic B and D 
decays |jTO|, radiative B — > K*'^ decay formfactor |ll| and B* —>■ Bn, D* —>■ Dn transition 
formfactors are calculated. The results of all works that have been devoted to the 
calculation of the formfactors confirm that, the predictions of the pole model are compatible 
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with that of the QCD sum rules method. Note that, in all these works the transition 
formfactors between the ground state mesons are considered. 

In connection with these observations, there follows an immediate question to whether 
the same agreement holds for the prediction of the formfactors resulting from the decays 
of the excited meson states. This article is devoted to find an answer to the question 
cited above, more precisely, we calculate the formfactors for the excited 1^ P wave decays 
Pi — s> P*7r, where Pi{Bi, Di) is the excited 1+ P wave meson state, and P* is the 1~ state. 
This article is organized as follows: In sect. 2, we derive the sum rules for the Pi P*7t 
transition formfactors. Sect. 3 is devoted to the numerical analysis and discussions. 

2 Sum Rules for the Pi P*7t transition formfactors 

For calculation of the sum rules formfactors for the Pi P*7t transition, we consider the 
following correlator, 

n^.fe,g) = ^ / d^xe'"'^ {n-{q)\d{x)^^Q{x)Q{0)^,^,u{0)\0) , (1) 

where, d7^Q and Q'ju'y^u are the interpolating currents of 1~ and l"*", respectively, Q is 
the heavy quark, q2 is the momentum of the 1~ meson and q is the pion's momentum. 
When the pion is on the mass shell, i.e., = m^, the correlator depends on two invariant 
variables g| and (g2 + g)^- In further calculations, we will set = 0. 

We start by considering the physical part of eq.(p. Before calculating the physical part 
of the eq.(l), few words are in order. According to the value of angular momentum of the 



light degrees of freedom (sf = |^), the heavy quark effective theory ||T2|, |T3|], predicts 



the existence of two multiplets, each including the l"*" meson: the first one contains 1^) 
and the second (2+, 1+) mesons. In terms of the conventional (2*+^)Pj states, the 1^ states 



defined above (physical states), are given by the following linear combinations |]T2|, |T4 



ii|> ^ \/f r^) /ipA) 

up ^ /ll'fl) - /fl'P.) . 

Therefore, the physical part of the eq.(l) must contain the contributions from both of the 
l"*" mesons, and it can be written as 

n _ V- {^~iQ)\dl,Q\lt){lt\Qlu7,u\0) 
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In |T3| , |TB|, it was shown that in the heavy quark hmit, tuq ^ rUg, 

(lt|Q7.75M|0) = . 

2 

The full theory can only bring small corrections to this result. For this reason the contri- 
bution of the it meson, in eq.(2), can safely be neglected. Therefore, in future analysis we 

2 

will consider the contribution of the it meson only, and for simplicity we denote it as Pi. 

2 

The matrix elements entering in eq.(0) are defined in the standard manner: 

{Pi\Ql,y75u\0) = -ifp,mp^e^{p) , 
(7r|J7,Q|Pi) = Foe,(p) + Fi(ge)(p + g), + F2(ge)(p-g), , (3) 

where /p^, eu and mp^ are the leptonic decay constant, vector polarization and mass of 
the excited 1^ meson, respectively. Fq, Fi and F2 are the formfactors that describe the 
Pi —>■ 71 transition. The contributions of the coefficients of Fi and F2 to the decay rate are 
proportional to 

ml + (mp^ - mp*)^ . 

From this expression it is obvious that, their contributions are small in comparison to the 
Fo's contribution. Our numerical analysis shows that the contributions of Fi and F2 to the 
decay width constitute about 20% of Fq, and hence, we shall neglect such terms. Thus, 



(^2 + qy - rrip^ \ rrip^ J 

and in further analysis we will pay our attention to the structure g^i, and denote its coeffi- 
cient by n. The physical part of the above relation takes the following form: 

^ ^ mf,^- (ga + qr -^^o s - {q2 + qf 

The first term on the right hand side is the pole term due to the ground state in the heavy 
channel, while the continuum and the higher order states are taken into account by the 
dispersion integral above the threshold Sq. 

Let us now consider the theoretical part of the correlator (|l|). This correlator function 
can be calculated in QCD at large Euclidean momenta (g2 + qY- In this case the virtuality 
of the heavy quark is of the order itlq — (g2 + qY-, and this is a large quantity. Therefore, we 
can expand the heavy quark propagator in powers of the slowly varying fields residing in 
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the pion, that acts as the external field on the propagating heavy quark. The expansion in 
powers of the external fields is also the expansion of the propagator in powers of a deviation 
from the light cone at = 0. The leading contribution is obtained by using the free heavy 
quark propagator in eq.(|l]) 



(2^ 



-ikx ' 



+ mQ) 



rrir 



Hence, 



d^xd^k e'('?2-*=)^ 



(7r(g)|g(x)7^(/{f + mQ)7j,75M(0)|0) . 



(5) 



We point out that, in eq.(^ and throughout the course of the present analysis, the path- 
ordered factor P exp {igs x^A^^ (ux) duj has been omitted, since in the Fock-Schwinger 
gauge x^A^ = 0, this factor is trivial. It follows from eq.(^) that, the answer is expressed 
via the pion matrix element of the gauge invariant, nonlocal operator with a light cone 
separation ^ 0. Following [jl^, the two particle pion wave functions are defined as: 

(7r(g)|(i(x)7^75M(0)|0) = -iq^,U f due 



tqux 



qx 



(finW + X gi[u] 

' du e''?""^2(^/) 



+ 







(7r(g)|J(x)75M(0)|0) 



du e 



iqux 



(6) 
(7) 



rriu + rrid Jo 

where (P-k{u) is the leading twist t = 2, ipp{u) is t = 3 and gi{u), (72 (^) are the r = 4 two 
particle pion wave function. Using eq.@ and eq.(0) after integrating over x and k, we get 
for the structure g^^: 

ml . (pp (u) , /•! ^,^ g2 (u) 



n(g2,g) 



—I 



rriu + TUd 



du- 



— h 2mQ I du—- 

A ^ Jo A 



-q2q 



du 



{gi {u) + G2 {u)) 

A2 



1 + 



2rr^ 
A 



"du^ 



where. 



A 

Goiu) 



^"q - (92 + quY 



g2{v)dv 



Since the higher order twist contributions are taken into account in eq.(^, the terms 
up to r = 4, which comes from the propagator expansion, must also be considered. The 
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complete expansion of the propagator is given in |T^, and it contains the contributions from 
the nonlocal operators qGq, qGGq , and qqqq. The only operator we consider here is qGq, 
because the contributions from the other two are very small, and hence they are neglected 
(for more detail see |18|). Under these approximations, the heavy quark propagator is 



given by the following expressions [y, [17 



Soix) 



1 



-ikx 



(2vr)4 



du 



rrir 



G^"'{ux)a^, + 



mg — k'^ 

Substituting second term in eq.(|^) into eq.(|l]), we get 

d'^x d^k du e'^'i^-''^ 



(2vr)^ (ml 



k^ 



(7r(g)|d(x)7^ uXpG^^{ux)-fx + 



{ml - k 
Using the following identities. 



1 

2 2 



"''^ ' GP\ux)a,x 7.75^(0)10) 



7m7a7i^ 



i {Qi^pIx - gf,xlp) + e^pAr7^75 , 



and the definition of the three particle pion wave functions |T7 



{n{q)\d{x)-fp-f5gsGai3{ux)u{0)\0) 



Ql3 ga^L 



Xaqfi 

qx ^ 



+ — {qaX/3 - qpXa) / Vanp\\{ai) e 
qx J 



qx y 

iqx{ai+ua3 



+ 



(9) 



(10) 



TT 



{q)\d{x)-ff,gsGap{ux)u{0)\0) 



X. 



aq^L 



qx 



+^{qaXf3 - qpXa) J Vai(p\\{ai)e 



qx / 



iqx{ai+ua^) 



+ 



we obtain from eq. (p!0|) , for the structure g^y 



n = iq2q I du I Va 



1 - 2M)v9||(ai) + ip\\{ai) 



2 ' 



(12) 



(13) 



^^a/SfiuG'^^ , is the dual tensor to 



where, Vai = daida2da3 5(1 — ai — 02 — 03), and Ga/s 
G'^'^, and y^'s are all twist four fuctions. Collecting all the terms (eqs.(|D and (||))for the 
theoretical part of the of the correlator function (HI), we get. 



n 



theor 



rriu + TUd 



A ^ 



A2 



du^-Af'dJ-^l^^^^ 
A Jo A2 

'l-2u)(f\\{ai) + (f\\{ai)\ 



1 + 



2m2 



A 



du / Vai 



where, 



A? 



Ai = niQ - [q2 + q {ai + uas)^ 



(14) 



The sum rule for the formfactor -fo(?2) obtained by equating the expressions 
for the invariant amplitude 11 (gg, {(I2 + ?)^)) which results as: 

Yjtheor 



and (III) 



fp^mp^Foiqi) ^ foo p^{qls)ds 



m 



A - (g2 + qf 



so 



s-{q2 + qY 



Invoking the duality prescripton and performing the Borel transformation in the variable 

iwii 



(g2 + qYi we get the following sum rules for the form factor -Fo(g| 



^0(92') 



fPimp, 

2 Jo 

xexp 



1 du 
— exp 



u 



M2 



ml - qlu 



^2[u,M\q' 



+ ^ / du 



m 



Pi 



(«! + -uas) 
ml - qlu 



rO(a;i + ua^ — 5) x 



M2 



mM2 



$3(M,M2,g2 



(15) 



where, 



u 



- ql 
So - qi 

1 — M , 



and. 



$9 



mu + md 
9i{u) 



2 2 
mi) — q^ 

mQLpp{u) H ^ (^^{u) + 



+2 
+2 



u 

G2{U) 

u 



1 + 



- ql 

uM^ 

- ql 



u 

^Qi^Q - qi) 



u^M^ 



+ 



(1 - 2u)Lp\\{ai) + (^||(a, 



- ql 



(«! + uas) M2 
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3 Numerical analysis and discussion 



Now we turn to the numerical calculations. In expression ([ISD , we use the following set of 
parameters: m^, = 5.738 GeV, mo, = 2.4 GeV, {so)b = 35 GeV^, {so)d = 8 GeV^. For 
leptonic decay constants fs^ and foi, we use the results given in |jl9|: fsi — 0.2 ±0.02 GeV 



and — 0.3 ± 0.03 GeV. The highest value of the Borel Parameter is fixed by 
imposing the condition that the continuum contribution is 30% of the resonance. With 
the help of this restriction, we calculate the maximum value of the Borel parameter for 
B{D) to be M^^^ ^ 20 GeV'^^c:^ 8 GeV"^). The minimal value of the same parameter is 
usually fixed by the condition that the terms proportional to the higher powers of -gj are 
negligible and it is found to have the value M^j„ ~ 8 G'eV^(:i^ 2 GeV"^) for B{D) mesons. 
Our calculations show that, the variation of within the above-mentioned limits, changes 
the result by less than 10%, which means that the dependence of the formfactor -fo(Q'2) on 
the Borel paremeter is quite weak. Note that, a similiar situation exists for the B ^ tt 
and D ^ IT transitions too (see for example [0,^1). Because of that, in our analysis we take 
= 15 GeV'^ for the Bi tx and = 4 GeV"^ for the Di tt transition formfactors. 
The maximum momentum transfer at which the sum rule (|15[) is applicable, is about 
15 — 20 GeV^ for B meson and 1 GeV^ for the D meson cases, respectively (for more detail 
see [0, 01). The explicit form of the pion wave function, used in eq.([T^) can be found in . 



The momentum transfer g| dependence of the formfactors F^{q2) and -^0^(52) are plotted 
in Fig.l and Fig.2. From these figures we observe that, 

^(f(92=0) = 1.2, 

F,^{ql = 0) = 1.5 . (16) 

Let us turn our attention to the pole model prediction for Fq^ and F^^ formfactors. In 
0], it is shown that the coupling constants Qe'e-k and gD*D-K fix the normalization of the 
formfactors of the B —>■ tt and D —>■ n transitions, respectively, within the context of pole 
model (see also, @, P). Using the pole model description in a similar manner, the formfactor 
^0(^2) can be expressed via the gp^p*-n coupling constant, that is calculated using the same 
correlator function (|l|) in |jT9|. Indeed, the formfactor F^^q^) defined by the matrix element. 



(7r(g)|J7^g|Pi(p)) =Fo(g2')eM(p) + ^i(92)(eg)(p + g)M + ^2(g2')(eg)(p-g)^ , (17) 
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is predicted to be (for the structure g^^) 

(i - ^ 

In deriving eg. (p!8D , we used the following definition, 

{nP*\P,)=gp^p,^{ee*) , (19) 

where e and e* are the 4-polarization vectors of the Pi and P* mesons, respectively. The 
dependence of the formfactors P(f (i?!) -Fo^((?2) I2 ^^e pole model (eq.(|l^)) are 
plotted in Fig.l and Fig.2, where we use gBiB*TT = 24 ± 3 GeV and gDiD*TT = 10 ± 2 GeV 
p!9| . From these figures we conclude that, in the overlap region both calculations agree. 



approximately, with each other. Quantitatively, at g| = 0, it follows from eq.(|T8D that 



F(f (g2 = 0) = 0.72 , 

F,''{ql = 0) = 1.20 . (20) 



If we compare eqs.(|T6D and (|20|), we see that, in the region for which uiq — g| > (1 GeV"^), 
{Q = b, c), the numerical agreement between the two approaches is better than 20% for D, 
but only within 35% for the B meson case. 

Finally we would like to point out that, these two approaches lead to absolutely different 
asymptotic behaviours of the relevant formfactors. Using the HQFT results, we get 

/p, ^/m = fp^ , fp* = fp* 

and, 

2m 

gPiP*n ^ -r- (see also [3J) . 



From eq.(|T8D it follows that, at mb — >• cxd we get, 

Fcf^^^ ~ {mp^n^r^' . (21) 

Then from eq.(p!5|) we see that, in this limit 

F^'^''^ ~ (mB(D))~'/' • (see also [20]) (22) 

In conclusion, we calculated the formfactors of the excited state l"*" mesons, namely, 
Bi —>■ 71 and Di ^ n transitions, in the framework of the light cone QCD sum rules 
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and compared our results with the pole dominance model predictions. The comparision 
elaborates that, the agreement between the two descriptions is rather good. This justifica- 
tion demonstrates that, the pole dominance model works quite good for the p-wave meson 
transitions, as it does for the s-wave ones. 
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Figure 1: The dependence of the formfactor Fo(g|), for the Bi ^ n transition, on gf|. 
The dotted Une corresponds to the hght cone QCD sum rules prediction and the sohd hne 
describes the pole dominance model prediction. 
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3 



2 - 




Figure 2: Same as in Fig.(l), but for the Di — > tt transition. 
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Figure Captions 

1. The dependence of the formfactor -Folo'i)) for the Bi ^ n transition, on g|. The dotted 
Une corresponds to the hght cone QCD sum rules prediction and the sohd hne describes 
the pole dominance model prediction. 

2. Same as in Fig.(l), but for the Di — > tt transition. 
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